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Quantum breathers are studied numerically in several electron-phonon coupled finite chain systems, 
in which the coupling results in intrinsic nonhnearity but with varying degrees of nonadiabaticity. 
As for quantum nonlinear lattice systems, we find that quantum breathers can exist as eigenstates of 
the system Hamiltonians. Optical responses are calculated as signatures of these coherent nonlinear 
excitations. We propose a new type of quantum nonlinear excitation, breather-excitons, which are 
bound states of breathers and excitons, whose properties are clarified with a minimal exciton-phonon 
model. 
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Much progress has been made toward understanding 
the physical consequences of nonhnearity over the last 
decade. In particular, recent developments concerning 
"breathers" 0| or "intrinsic local modes" (ILM) , sug- 
gest that energy focusing is prevalent in both classical 
1^-^ and quantum [^[jsj discrete nonlinear nonintegrable 
frameworks. The existence and stability of multiquanta 
bound states (breathers, ILMs) are now established for a 
wide variety of discrete classical models with prescribed 
nonhnearity. Recently, we demonstrated that this prop- 
erty also persists in a discrete quantum nonlinear lattice, 
and exhibited some distinctive observable signatures in 
terms of spatial-temporal correlations However, of 
more profound concern is the typical origin of effective 
nonhnearity in quantum systems, namely, through the 
coupling of two or more fields. Adiabatic slaving of fields 
usually results in nonlinear Schrodinger models. Realis- 
tically, however, nonadiabatic effects must be considered 
and the influences of nonhnearity and nonadiabaticity 
are inevitably interrelated. Here we consider examples of 
electron-phonon (e-ph) coupled models frequently used 
to describe organic and inorganic correlated electronic 
materials, for which an adiabatic treatment of breathers 
is inadequate for various physical observables |^ . The in- 
teractions of electrons with the lattice and among them- 
selves provide sources of nonhnearity and strongly influ- 
ence electronic, optical and structural properties. Here, 
we use numerical approaches restricted to finite chains 
but without any adiabatic approximation. We find that 
quantum breathers of two types can exist in the e-ph sys- 
tems considered — those near electronic ground states, 
and photoexcited breathers (termed "breather-excitons" 
below) . 

We consider a Holstein-Hubbard (HH) tight-binding 
model Hamiltonian [p^ of an e-ph coupled system: 

= ^ -to(cLc,+ l,. + H.C.) (1) 
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Here, cj^. (cio-) and bl (bi) are the creation (annihilation) 
operators of electrons and phonons, respectively, to, U 
and V are the electron kinetic energy, and the on-site 
and nearest-neighbour Coulomb repulsions, loq is the 
bare phonon frequency. The form of the e-ph coupling (A 
term) is that used in the Holstein (H) model iQ . We use 
models with one electron per site {i.e. a ^-filled band). 

Our numerical approach mainly consists of exact di- 
agonalization of Hamiltonian matrices represented in 
Hilbert spaces defined by appropriately selected basis 
functions, and the analysis of characteristics of quan- 
tum breathers via various dynamic correlation functions. 
We deal fully with nonadiabaticity and e-e correlations; 
the only approximation is the truncation of the infinite 
phonon Hilbert space. Because wavefunction information 
is needed to identify a quantum breather state which is 
not necessarily low-lying in the complete eigenspectrum, 
efficient exact diagonalization techniques have been de- 
veloped to handle large-scale matrices; these are also de- 
signed to suit parallel computer architectures In 
Fig. |l|, we show a typical total density-of-states (DOS) 
of a HH model, calculated using the kernel polynomial 
method ||l^ which has a finite energy resolution. The 
regions of interest are indicated (see text below Figs. ^ 
and ^). Comparing to the Hubbard model DOS, we real- 
ize that the denser the states the more difficult they are 
to numerically distinguish, and the higher the numeri- 
cal efficiency required. Also, since our model systems 
are finite, it is not easy to eliminate finite-size effects 
jisf . However, our main results remain valid because the 
breather excitations of interest here are intrinsically of fi- 
nite size. There are also many realistic material contexts 
which are of finite-size (e.g. conjugated oligomers). 

The eigen-energies of the system (|l|) are shown in Figs. 
H and|| for 6-site Holstein, and HH models, respectively. 
The results in Fig. |^ were obtained using a basis set 
of Debye phonons and Bloch electron functions, while 
Fig. □ results use Einstein phonons and Wannier func- 
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tions. The degrees of softening of the k = tt mode (Kohn 
anomaly) indicate that the systems in Figs. 2(a) and 
3 are weakly e-ph coupled, whereas that in Fig. ||(b) is 
more strongly coupled. The phonon dispersion is caused 
by effective nonlinear ph-ph coupling via the e-ph cou- 
pling. Because of discrete translational invariance, direct 
examination of the wavefunction amplitudes and single 
phonon operator expectation values do not reveal the 
presence of quantum breathers |8|. Rather, we examine 
various static and dynamic correlation functions of lattice 
displacements, including: Uj{t)=J2i{k\ Ui{0)uj{t) |fc), 
where Ui is the phonon displacement operator. The static 
correlation functions (t=0) probe the spatial localization 
in a given eigenstate, while the dynamic counterparts 
probe the temporal coherence. These two properties are 
the distinguishing characteristics of quantum breathers 
10,1). The correlation functions Uj{t), as in Fig. ^, show 
that there indeed exists a band of particle-like states in 
each of the above examples, and that these states pos- 
sess short but finite spatial and temporal correlations, 
whereas all other states are extended [|l^ . A typical case 
is shown in Fig. ^(b), where the particle- like band pos- 
sesses a large binding energy separating it from the con- 
tinuum bands above, and the anharmonicity and corre- 
sponding localization are strong (the correlation length 
is approximately 3 lattice constants, see Fig. ^ pof . 

The existence and properties of quantum breathers de- 
pend not only on the effective nonlinear ph-ph coupling 
(as in ||]) but also on the nonadiabaticity (here, con- 
trolled by the ratio Tiujo/to). We find that the stronger 
the e-ph coupling and the adiabaticity, the more eas- 
ily quantum breathers form p^ . Furthermore, these 
breather states can survive strong e-e correlations (Fig. 

, although we find that in the lowest part of the ex- 
citation spectrum strong e-e correlations tend to induce 
extended magnetic excitations jl^]. We will see below, 
however, that strong e-e correlations do not necessarily 
destroy short-correlation length breathers in other spec- 
tral regions. 

To identify physical consequences of the breather ex- 
citations, we study the first-order optical response func- 
tion which can be measured in one-photon experiments. 
Fig. ^ shows the zero-temperature infrared and electronic 
optical absorptions of a HH system. Within this nona- 
diabatic approach, we account for all the electronic and 
phonon polarizabilities within one photon perturbations. 
This is a step toward understanding accumulating exper- 
imental results in ultra-fast time-resolved (nonadiabatic) 
and nonequilibrium measurements jl^ . First, we observe 
that in the infrared region, several prominent peaks di- 
rectly indicate the anharmonicity, including multiphonon 
side-bands. Comparing to Fig. |^, the contribution to 
the infrared absorption from the breather states (marked 
"&") has comparable intensity to that from surround- 
ing extended states. Second, below the main absorption 
peak (marked "exciion"), there exist a series of spectral 



features which cannot be explained in an adiabatic de- 
scription. To understand their origin, we added small 
amounts of disorder in the electron or phonon degrees of 
freedom {e.g., the parameters U or Titoa). We found that 
only phonon disorder changes the spectral features near 
the "excitori^ edge p^], indicating their phonon origin. 
Third, in addition to the rich structure in the region of 
electronic absorption, we observe a series of bands be- 
tween the "ea;ci<on" edge and the continuum" bands. 
Among them, we notice one band ("&-e") whose posi- 
tion and intensity relative to the other phonon sidebands 
does not change under the perturbation of phonon disor- 
der jl^. This indicates that this particular band is due 
to some a relatively stable excited configurations which 
are more localized. This evidence suggests a possible 
new type of bound state which is higher in energy than 
the lowest excitons and more localized than the nearby 
states — we will term these breather- excitons. Below we 
show that they possess some of the characteristics of the 
"ground state" quantum breathers identified in the low- 
lying part of the eigenspectrum [Fig. H). 

Physical intuition suggests the possibility of a bound 
state of excitons and breathers, i.e. a hot (dressed) exci- 
ton or a photoexcited breather. Excitons exist primarily 
as electron-hole pairs bound by both the Coulomb in- 
teractions and e-ph coupling. The latter factor slows 
down the exciton motion and tends to dress the exciton 
with phonons and breathers, e-e correlations provide an- 
other energy region from which the breather can be ex- 
cited, and contribute additional nonlinearity enhancing 
the breather formation. Furthermore, with the excited 
oscillating dipole moments inside the breather-exciton 
bound states, they will strongly absorb photons, as an 
electronic exciton does. 

To support the above arguments [p^ , we introduce a 
minimal exciton-phonon model, describing an electron 
and hole interacting with each other and with phonons: 

H = Y, -te {^el e,+i + H.c. ) + 

i i 

hi h,+i + H.c. ) -Y^f" 

i i 

+ E V{i-j)ele,hlh, 

i i 

+ Y Ulh + l) , (2) 

with V{i—j) — —U (as i = j); ~V (as \i—j\ — 1); 
and (otherwise). e| (e^) and /i| (hi) are creation 
(annihilation) operators for the electron and hole, re- 
spectively. ee(h), Vii - j) and Ae(h) are the elec- 
tron (hole) hopping integrals, on-site energies, electron- 
hole attractions and e(/i)-ph couplings strengths, respec- 
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tively . Fig. ^ shows the low-lying part of the eigen- 
spectrum, corresponding to the "eiciton" region in Fig. 
||. The lowest band corresponds to the lowest exciton 
states in Fig. |^ and shows dispersion. They are exci- 
tons dressed by a static lattice distortion of condensed 
phonons. The breather-exciton states are also strongly 
localized, but dressed with a dynamic phonon wavepacket 
(weak dispersion) and exhibiting more coherent inter- 
nal motions. (We speculate that the excitonic internal 
frequency resonantly traps a breather.) The remain- 
ing states are e-h pairs coupled to extended anharmonic 
phonons. Breather-excitons are more polarized than the 
ground state breathers, leading to the stronger optical 
absorption seen in Fig. ^. 

In summary, this study is an important extension of 
previous investigation of quantum breathers in nonlinear 
lattices 1^ . It demonstrates their existence in more gen- 
eral e-ph systems where nonlinearity is self-consistently 
generated by the e-ph (and e-e) coupling rather than 
being inserted by hand. Thus the combined influences 
of nonadiabaticity and nonlinearity are incorporated. 
Clear features from breathers are found in optical proper- 
ties. The concept of breather-exciton states was proposed 
and demonstrated for a simplified exciton-phonon model 
without any adiabatic approximations. Investigations are 
underway to correlate our findings with nonlinear op- 
tics and time-resolved spectroscopy, including the com- 
peting timescales associated with exciton formation, ex- 
citonic self-trapping, and breather-exciton formation, as 
functions of nonlinearity and nonadiabaticity [ p^ . 
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FIG. 1. Total DOS of a 6-site Holstein-Hubbard chain 
with to = 1.00, uj = 0.40, A = 0.31, U = 8.10, V = 1.11, 
Mph = 4, and X = 1 638 400. (Mph and M stand for phonon 
truncation and the total dimension of the matrix.) The inset 
shows the DOS over the full energy range, with the corre- 
sponding pure Hubbard model (gray line) results for compar- 
ison. 



FIG. 2. Eigenspectrum (low-lying sector) of a 6-site Hol- 
stein chain with fo = l-00, 1.^=0.4, U^V^O.OO and (a) A = 0.31, 
Mph(fc = ±f , ±^) = 5, Mph(fc = n)=6, Mph(fe = 0) = 1, 
M = 1500000; (b) A = 0.61, Mph(A; = ±f , ±^) = 6, 
Mph(fc = tt) = 10, Mpi,{k = 0) = 1, M = 5184000. The 
breather states are labeled with solid triangles, and the re- 
maining ones with open triangles. The ground state energy 
has been subtracted from E^. 
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FIG. 3. Eigenspectrum of the 6-site Holstein-Hubbard 
chain in Fig. |^. The notation is as in Fig. ^. The ground 
state energy is 4.14to 

FIG. 4. The spatial correlation function Uj{0) of the sys- 
tem in Fig^(b). The solid line is for a breather state (labeled 
B in Fig. 0(b)), while the gray line is for an extended state 
(E). 



FIG. 5. The zero-temperature infrared and electronic op- 
tical absorptions of the 6-site Holstein-Hubbard system in 
Fig. ^ (solid line) and the related 6-site Hubbard model with 
to = LOO, U = 8.10, V = 1.11 (gray line). AU spectra are 
broadened by a Lorenzian of width 0.005. 
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FIG. 6. Eigenspectrum (low-lying sector) of the 6-site 
minimal exciton-phonon model (see text) with te = 1.00, 
th = 0.67, ee = 4.00, eh = -4.00, loq = 0.15, U = 18.0, 
V = 2.75, Ae = 0.87, = -0.27. Mph = 6, X = 1679 616. 
The breather-exciton states are labeled with solid triangles, 
the dressed excitons with diamonds, and the remaining ones 
with open triangles. The ground state energy has been sub- 
tracted from Ek- 
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Figures in "Quantum Breathers in Electron-phonon Systems" by 
W. Z. Wang et al, November 1997. 




FIG. 1. Total DOS of a 6-site Holstein-Hubbard chain with to = 1.00, ui = 0.40, A = 0.31, U = 8.10, V = 1.11, Mph = 4, 
and M = 1638 400. (Mph and M stand for phonon truncation and the total dimension of the matrix.) The inset shows the 
DOS over the full energy range, with the corresponding pure Hubbard model (gray line) results for comparison 
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Figures in "Quantum Breathers in Electron-phonon Systems" by 
W. Z. Wang et al, November 1997. 
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Figures in "Quantum Breathers in Electron-phonon Systems" by 
W. Z. Wang et al, November 1997. 
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FIG. 3. Eigenspectmm of the 6-site Holstein-Hubbard chain in Fig. 1. The notation is as in Fig. 2. The ground state energy 
is 4.14to 
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Figures in "Quantum Breathers in Electron-phonon Systems" by 
W. Z. Wang et al, November 1997. 
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FIG. 4. The spatial correlation function 11^(0) of the system in Fig. 2(b). The solid line is for a breather state (labeled B in 
Fig. 2(b)), while the gray line is for an extended state (E). 
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Figures in "Quantum Breathers in Electron-phonon Systems" by 
W. Z. Wang et al, November 1997. 
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FIG. 5. The zero-temperature infrared and electronic optical absorptions of the 6-site Holstein-Hubbard system in Fig. 3 
(solid line) and the related 6-site Hubbard model with to = 1.00, U = 8.10, V = 1.11 (gray line). All spectra are broadened by 
a Lorenzian of width 0.005. 
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Figures in "Quantum Breathers in Electron-phonon Systems" by 
W. Z. Wang et al, November 1997. 
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FIG. 6. Eigenspectmm (low-lying sector) of the 6-site minimal exciton-plionon model (see text) with te = 1.00, tu = 0.67, 
ee = 4.00, eh = -4.00, loo = 0.15, U = 18.0, V = 2.75, Ae = 0.87, en = -0.27. Mph = 6, M = 1679 616. The breather-exciton 
states are labeled with solid triangles, the dressed excitons with diamonds, and the remaining ones with open triangles. The 
ground state energy has been subtracted from Ek- 
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